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ABSTRACT 

The geomagnetic electrokinetograph (GEK) was used to measure 
surface currents near the center of Monterey Bay during six separate 
24-hour periods from May through July, 1972. An average of 244 current 
vectors were derived for each cruise. The mean currents from these 
cruises are all southerly and ranged from 4.1 cm/sec to 2 0.4 cm/sec. 

The average of these mean currents is 12.1 cm/sec toward 163°T. These 
values were compared with individual currents derived from dynamic 
topographies from the same period. 

Diurnal and semi-diurnal variations of the current were studied 
after subjecting the data to a fourier analysis. It was concluded that 
there must be at least an indirect coupling of the ocean currents with 
the semi-diurnal tide at the data point. The diurnal component also is 
important; it may be tidal or inertial, or merely related to the passage 
of the sun. 
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I. INTRODUCTION 



I. INTRODUCTION 

This study was designed to use the Geomagnetic Electrokinetograph 
(GEK) to investigate the time variability of surface currents at one point 
in Monterey Bay. 

Recent studies by Smith (1971) and McKay (1970) utilized the GEK 
to measure surface currents, but these measurements were taken singly, 
over large areas, and at different times. In order to determine the nature 
of the surface currents in Monterey Bay it was necessary to get data 
which was independent of spacial changes and of a sufficiently long 
duration to show tidal and inertial effects. The technique followed was 
to continually run squares with the GEK for approximately 24 hours. The 
squares averaged 0.5 nautical miles on a side and were repeated about 
every 24 minutes. Electrode zeros and current components were com- 
puted cyclically around the square. 

Six experiments were conducted between 10 May and 22 July, 1972 
at a point near the center of Monterey Bay (36 o 46'N/l22 o 01’W) in water 
440 fathoms deep (Figure 1) . The Naval Postgraduate School (NPS) 
research vessel ACANIA was used to tow the GEK array. Analysis of 
the data was facilitated by the use of the NPS IBM 360 computer. Causes 
of variability of surface currents were sought in tidal, wind, inertial, 
and seasonal effects . 
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II. THE ELECTROMAGNETIC METHOD OF 
MEASURING OCEAN CURRENTS 



A. DEVELOPMENT 

When a part of the ocean moves relative to the earth's magnetic 
field, a voltage is induced in the sea water. Under certain conditions 
a voltage equal to the induced voltage may be measured between two 
identical electrodes introduced into the sea. Under other conditions, 
part or all of the induced voltage is dissipated through resistance in the 
sea water. 

Faraday is credited with the first mention of and experimentation 
with this phenomenon. His experiments failed; however, in the late 
1800's broken trans -Atlantic cables showed the potential differences 
he had sought and led to the early 20th century experimentation in 
Dartmouth Harbor by Young, Gerrard, and Jevens (1920). 

The electromagnetic method of towed electrodes is a continuation 
of this theory and is a result of work done between 1946 and 1954 by 
Longuet-Higgins and Barber (1946), Stommel (1948), von Arx (1950), 
and Longuet-Higgins, Stern, and Stommel (1954). The technique owes 
its name. Geomagnetic Electrokinetograph (GEK) , to von Arx. Using the 
GEK method, surface currents can readily be measured from a ship underway. 

B. GEOMAGNETIC ELECTROKINETOGRAPH (GEK) 

1 . Theory of Operation 

The GEK is made up of a recording potentiometer, a two- 

conductor, neutrally buoyant cable, and a matched pair of silver-silver 
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chloride electrodes of the order of 100 meters apart. The potentiometer 
records the difference in potential between two electrodes. A long cable 
is required to allow the electrodes to be towed beyond the magnetic 
effects of the ship's steel hull. 

An extensive discussion of electromagnetic theory and the ocean's 
movement through the earth's magnetic field may be found in von Arx 
(1950) and Longuet-Higgins , Stern, and Stommel (1954) (hereafter referred 
to as L-HSS) . The following summary of L-HSS's treatment of this sub- 
ject is by McKay (1970) with large portions of the development based on 
work by Curtin (19 70). 

The electrical field intensity along a conducting path that moves 
through a magnetic field is the difference between the two potential 
gradients. One potential gradient can be shown by application of 
Faraday's Law of Induction to be V x H as is done by L-HSS. In this 
expression, V is the vector velocity with which the conducting path is 
moved and H is the magnetic field intensity vector. The other potential 
gradient is given by Ohm's Law as />_i , where /) is the resistivity 



and_i the electric current density vector along the conductor. If the 
potential is denoted ^3 , the negative of the electric force field in 

the moving conducting path can be expressed as 



When the electrodes are deployed, two distinct conducting paths 
are present: (a) the sea water, and (b) the long conductors to which the 
electrodes are attached. The potential in each of these paths is the 




( 1 ) 
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difference of the two terms of equation (1) . In the following discussion, 
a right-hand co-ordinate system is utilized with the following conven- 
tions: the positive x-direction is that of the ship's head; y is the 
horizontal co-ordinate perpendicular to the course steered and z is 
positive up. 

Consider first the sea water path. When the electrodes are towed 
directly astern of the ship, as shown in Figure 2, they sense the poten- 
tial difference along a path parallel to the direction of the ship's head. 
The component of potential gradient generated in this direction due to 

the (VxH) term is (V H - V H ) , where V is the component of water 

y z z y y 

velocity perpendicular to the course steered (set of ship), V is the 

z 

vertical component of water velocity, H is the vertical intensity of the 

z 

earth's magnetic field, and H is the horizontal intensity of the geo- 
magnetic field in a direction perpendicular to the ship's course. If V 

z 

is assumed negligible in comparison with V , the measurable potential 



gradient generated by the water motion is V H ; that is , the potential 

Y z 

gradient contributed by the (VxH) term is due to the component of water 
velocity perpendicular to the ship's course. This water motion causes 
the ship's set (and the set of the electrodes) and is independent of the 
ship's velocity in the direction of the interelectrode line. The second 
potential gradient in the sea water path, yO _i, is a direct function 
of the electric current density along the path between the electrodes. 

The electric current density along this path depends on the electrical 
characteristics (potential gradient’s resistances) of all other surrounding 
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paths, since it is only through the surrounding paths that a closed 



electric-current-carrying circuit can be completed . 

Consider next the part of the closed conducting path formed by the 
two wires and the voltmeter. It can be seen, by application of reasoning 
identical to that used for the sea water path, that the velocity -induced 
potential gradient due to the ( V x H) term in equation (1) is again V H . 



due to the low resistance of the wire and the small current drawn. 

The potential differences sensed by the voltmeter and continuously 
recorded as GEK data are those developed in the wires leading to the 
electrodes and in the sea water between the electrodes. Traversing this 
path, as shown in Figure 2, and adding the potentials yields the following: 



where E is the emf per unit distance and A, B, and C are the points shown 
in Figure 2. The sumbol AB represents the distance between A and B. 

The term in brackets represents the potential gradient in the sea water 
between the electrodes; the remaining terms represent the induced emf 
in the indicated sections of the wire . Equation (2) reduces to: 



where the direction is that of the electrode line, or, in the case of 
electrodes towed directly astern, the direction of the ship's head. 



y z 



The second potential gradient is essentia Uy zero in this path 




( 2 ) 



+ (-V H )BA 
y z 




(3) 
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riguro 2. GEX COMPONENT 
ARRANGEMENT 



For most currents in the deep ocean the potential gradient is 



completely short-circuited by the deep layer of relatively still water 
underlying the moving surface layer. In this case (^) = 0 and 

d = Vx H. Under these conditions, if the interelectrode distance 



is s units, the horizontal potential difference measured by the voltmeter 



is: 

E = (ksH z V h ) x l(f 8 (4) 

where V, is the horizontal component of water motion perpendicular to 
h 

the interelectrode line and k is a constant near unity, to be discussed 
in Section III C . 

The positive direction of E is taken as that of the ship's head 

and the current velocity is determined by adding the electric vectors, E, 

obtained on two courses, at right angles. The following example indicates 

the relationship between the direction of the current component measured 

and that of the interelectrode line: 

If V is positive toward the east, H is directed downward, 
h z 

and the interelectrode line is parallel with a meridian of 

longitude, E will have a positive polarity at the southern 

electrode and a negative polarity at the northern electrode. 

Since II is positive downward in the northern hemisphere, the 
z 

following convention must be observed to obtain the true direction of 
the current velocity: In order that the deflection recorded by the in- 
strument be in the same direction as the set of the electrodes, the 
conductor leading to the more distant electrode must be connected to 
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the voltmeter input terminal which, when made positive, gives a de- 
flection of the recorder current indicator to the right. The necessity of 
observing this convention can be seen from the example given above. 

This convention must be reversed in the southern magnetic hemisphere, 

where H is directed upward, 
z 

2 . Assumptions 

As has been stated above, the electric current density, _i, 
is a function of the conducting path in which it exists. Thus, with an 
assumption of uniform resistivity, the current density integrated along 
any line is directly dependent on the emf's in the available sea water 
paths surrounding that line. These emf's as indicated by Equation (1), 
are directly related to the water velocity in the region. Therefore, any 
electric current density measured at the surface is a reflection, in part, 
of the velocity field in the entire water column below the surface and in 
adjacent waters. 

Recognition of this fact is essential for a meaningful inter- 
pretation of GEK data. It is only in the case of true oceanic currents 
that GEK data, unsupplemented by information derived from other sources, 
can be considered representative of the true surface current velocity. In 
this case the following conditions generally prevail: 

(i) the vertical component of water velocity is negligible 
in comparison with the horizontal component, 

(ii) the flow is broad compared with the depth of the water, 

(iii) the current is steady in the x-direction, while a current 

"fix" is being obtained, 
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(iv) a deep layer of still water underlies the moving surface 
layer. 

If condition (iv) is not obtained, the voltage is usually less than 
the theoretical and it is customary to multiply by a correction factor k 
determined by a simultaneous measurement with a cun'ent meter. 
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III. EQUIPMENT 



A. ELECTRODES 

A complete operating assemblage of GEK equipment was inherited 
from Smith (1971) at the onset of this project. 

The matched set of silver-silver chloride electrodes was supplied 
by Woods Hole Oceanographic Institute and attached to the cable 100 
meters apart by direct soldering. Von Arx (1950) describes the process 
by which electrode pairs are matched. The electrodes were padded and 
housed in a water permeable tube as shown by Smith (1971). 

During the three-month period of data collection, it was not neces- 
sary to change the electrode installation. 

To eliminate the electrochemical potential, the "zero point" must 
be obtained and subtracted out of the actual current vectors. This is 
normally done by reversing the electrodes end for end. In this case, the 
ship used a square sailing plan and the zero point was taken from opposite 
sides of each square. This is described in detail in Section V A. 

B . CABLE 

A neutrally buoyant cable (Type B-3088) 329 meters long, made by 
the Vector Cable Company, was utilized during this investigation. The 
cable's function is to transmit to the recorder the electrode-generated 
signals from far enough astern (1.5 - 3.0 shiplengths) of the ship that 
the electrodes are not affected by the ship's magnetic field and wake. 
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C. RECORDING POTENTIOMETER 

The Geomagnetic Electrokinetograph Model V recording potentiometer 
based on a Speedomax Type G recorder was used to record the potential 
differences delivered from the electrodes by the cable. A photograph 
and a thorough description of the operation of the GEK Model V is con- 
tained in H. O. Publication 607. 

The vertical intensity of the earth's magnetic field in gauss for 
the ship's location was taken from H. O. Chart Number 1702 and the 
rounded value was set into the recorder to the nearest tenth. A multi- 
plication factor was later used to adjust the data to the exact magnetic 
vertical intensity at the data point. The multiplier includes compensation 
for annual variation, and is the ratio of the rounded vertical intensity to 
the correct local magnetic intensity. In this instance it was 0.91 
(0.40/0.44). 

As noted by Smith (1971) it is essential to isolate the power for 
the recorder from the ship's ground. This is done with an isolation 
transformer which must not be by-passed by the ground lead of a 3-wire 
power cord. Hence, a conventional two-prong plug is necessary. 

The recorder also requires a signal multiplier setting. This k 
value is the ratio of the water speed observed by a non-elcctromagnetic 
method to the water speed measured by the electromagnetic method. An 
open ocean value of 1.04 is suggested by von Arx (1950) with the caution 
that higher values might be needed in areas where tidal currents and 
water shear occur. 
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The k factor was set at 1 . 10 following the recommendation of 
H. O. Publication 607 (1970 Ch. 1) for water deeper than 75 fathoms. 
This factor was derived from experience in open ocean conditions along 
the Atlantic coast. An unusual structure like Monterey Canyon may 
behave anomalously, a possibility which will be considered in the dis- 
cussion of accuracy (Section V D) . An incorrect choice of k does not 
lead to an error in direction. 
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IV. SHIP PROCEDURES 



A. PLATFORM 

The research vessel ACANIA was acquired by the Naval Postgraduate 
School in 1971 and refitted as an oceanographic research vessel. She 
was launched in 1929 and for many years took part in research activities 
under the Stanford Research Foundation. 

R/V ACANIA' s length is 126 feet, maximum speed is ten knots, and 
she carries a crew of five. Of the standard equipment aboard ACANIA, 
the Decca RM 416 radar was the most important to this study as noted in 
Section IV C. The UQN - 1 depth finder was also used briefly to confirm 
charted depths at the data point. 

B. SAILING PLAN 

Both clockwise and counterclockwise squares were traversed during 
the cruise series, allowing the ship to remain within 1000 meters of the 
data point in calm weather and within 1500 meters during heavier weather. 

Ship's speed was held at five knots for all cruises, but the time 
on each leg was varied between four and nine minutes, depending upon 
wind and current conditions . In order to maintain position within 1000 
meters of the data point the downwind leg might take four minutes , while 
the upwind leg (opposite side of the square) might take eight minutes. 

All legs were at least four minutes long to insure that the cable had 
straightened out and a good signal was being received. 



22 



FIGURE 3 . SAILING 
FIRST VECTOR: 



SECOND VECTOR: 
THIRD VECTOR: 





V4 


' / 


V3 


DATA 


VI 




POINT 




4 


V2 





VS 



v 



PLAN AND COMPONENT COMPUTATION 

ZERO = VI + V3 
2 

E-W COMPONENT = Vl-ZERO 
N-S COMPONENT = V2-ZERO 

2 2 1/2 

CURRENT MAGNITUDE = [ (E-W) + (N-S) ] 

CURRENT DIRECTION = arcsin E-W 

MAGN 

ZERO = V2 + V4 
2 

E-W COMPONENT = V3-ZERO 
N-S COMPONENT = V 2 -ZERO 

ZERO = V3 + V5 
2 

E-W COMPONENT = V3-ZERO 
N-S COMPONENT = V4-ZERO 



\ 
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Only cardinal headings were steered during this series of cruises. 



The sailing plan is shown in Figure 3. 

C . NAVIGATION 

The data point for this project was located close to the center of 
Monterey Bay over the Canyon. Navigation was by radar at ranges of 
roughly 11 miles from well-defined landmarks. These points are noted 
on Figure 1. Visual bearings were also used to augment the radar range 
information. The positioning accuracy is estimated to be + 0 . 1 nm. 
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V. DATA COLLECTION, REDUCTION, AND PRESENTATION 



A. SURFACE CURRENTS 

1 . Data Collection 

During each cruise approximately twenty-four hours of con- 
secutive surface current readings were obtained. All data was taken at 
one point near the center of both Monterey Bay and the Monterey Canyon. 
An average of 244 current vectors were obtained from each of the six 
cruises . 

The GEK recording potentiometer was continuously manned 
in order that the strip chart could be annotated and the potentials logged 
for each leg of the sailing plan. The ship's speed was five knots with 
each leg of the square being between four and nine minutes in duration. 
The average duration was six minutes. Longer and shorter legs were 
run to compensate for the ship's drift. 

The primary items logged on the data sheet were: ship's 
heading (°T) , velocity (cm/sec), and time. Weather, position, depth, 
speed, and other pertinent data were also logged on each data sheet. 

The strip chart was marked and the time was noted for each 
change of course. The date, segment of cruise and cruise number were 
also recorded on the chart. 

Breaks in the data of about ten to fifteen minutes occurred 
when the direction sailed around the square was changed from clockwise 
to counterclockwise or vice versa. Several breaks in the data sets of 

up to forty-five minutes occurred due to steering casualties. 
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Data Reduction 



2 . 

The IBM 360/67 computer at the Naval Postgraduate School 
computer center was used to reduce the rough data to current magnitudes 
and directions . 

As can be seen from Figure 3 , there is an interdependence 
of current vectors in each square since one of the current vector com- 
ponents is calculated from a leg of the square which will take part in the 
next zero computation. This procedure allowed for the calculation of 
approximately ten current vectors per hour. 

The outputs of the data reduction computer program included 
east-west (E-W) and north-south (N-S) components, time, and the mag- 
nitude and direction of each current vector. This computer program is 
listed in Appendix A with suitable comments and explanations. Appendix 
B lists sample current data for Cruise I. 

Sharply alternating speeds were present in the data after 
reduction. These variations were removed by averaging alternate pairs 
of data. A full discussion of these variations is to be found in Section V 
D. 

3 . Data Presentation 

To produce a clear presentation of current magnitudes versus 
winds and tides, a computer program was written which utilized the 
plotting package for the CALCOMP Model 765 plotting equipment along 
with the IBM 360 computer. 
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For each cruise this graphing program produced the plots of E-W 
and N-S components and resultants versus time. Corresponding tide 
heights and wind speeds were plotted on the same graph. Northerly and 
easterly sets were taken to be positive. Both unaveraged and averaged 
data were handled in this fashion. The graphing program is included 
after the reduction program in Appendix A. 

The data is also presented in a progressive vector diagram format 
for each cruise. The progressive vector diagram again utilized the plotting 
package and may be found with explanations in Appendix A. 

B. WINDS 

Estimations of wind speed and direction were recorded aboard R/V 
ACANIA at various times (normally every two or three hours) during each 
cruise. These estimates were compared with the hourly readings logged 
at the U. S. Coast Guard lighthouse on Point Pinos. A composite of 
these two sources forms the basis for the wind magnitude and direction 
utilized for each cruise. Appendix C contains the resulting hourly wind 
data for Cruise I. 

C. TIDES 

Tide records from Pier Number 1 at Monterey were made available 
by the NPS Oceanographic Department. This data was used to plot 
directly the tidal curve for each cruise. These tides compare favorably 
with the high and low tidal heights as calculated from the U. S. Coast 
and Geodetic Survey Tide Tables. The hourly tidal heights utilized arc 



contained in Appendix D. 
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A tide prediction computer program was also written, based on 



harmonic constants for San Francisco and corrected to the data point in 
Monterey Bay. The results from this effort did not correlate well with 
either the tide gauge measurements or the tide tables, and were, there- 
fore, not considered further. 

D. DISCUSSION OF ACCURACY 

In discussing the accuracy of the current data the following poten- 
tial errors must be considered: 

(i) Errors in the observed voltage due either to reading in- 
accuracies or to the failure of the trace to steady to a 
constant value after a turn. 

(ii) Errors due to inadequate depth or effects of the canyon 
walls not compensated for by the k factor. 

(iii) Errors due to substantial velocity gradients across the 
square . 

Since the trace could be read to + 1 cm/sec, and at least four 
minutes were always given to steady out this trace, the errors in Group 
(i) are small. They are not systematic and in the long term average to 
zero . 

It is not known how great the errors of Group (ii) may be. There 
may not be a thick layer of static water near the bottom (current shears, 
internal waves may exist) and the lateral dimensions of the parcel of 
deep water are relatively small. Because of the sharply sloping walls 
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of the canyon, bottom effects may affect different legs unequally and 
cause aberrations. In these circumstances it is unwise to be categorical 
about the adequacy of using k = 1.10. On the other hand, the speeds 
measured are plausible and not notably at variance with the few available 
measurements obtained in various locations in Monterey Bay. 

The sharply alternating speeds described in Section V A 2 may best 
be explained by the existence of large velocity gradients across the 
square. These gradients may be caused by variations of depth or by 
actual currents or eddies. This is contrary to the assumption that the 
current is broad and therefore the same on both sides of the sailing plan. 
This breakdown, coupled with the circular and redundant method by which 
the currents were computed, produces oscillations in a sequence high- 
mcdium-low-medium-high in the N-S and E-W components of current 
and in the resultant. (See Appendix C for an example.) If the current 
component actually present on one N-S leg has a magnitude E and on 
the opposite side has the magnitude E + ZA E, and likewise the re- 
maining two sides of the square have magnitudes N and N + ZA N , 
then it may be shown that the peak-to-peak range of the resulting zig- 
zag will be the absolute value of A E - A N. (Sec Tigure 3 for 
method of computation.) A maximum peak-to-peak range of 60 cm/sec 
occurs in the E-W component of Cruise I, though most peak-to-peak 
variations are of the order of 30 cm/scc. These variations were mostly 
removed from the data by averaging alternate pairs of values. Variations 
after averaging are of the order of 5 cm/sec. The process of averaging 
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alternate pairs would produce a speed corresponding to some point within 



the square, if the first explanation of the alternation is valid. It can be 
seen from comparing the graphs of averaged and unaveraged data, in 
Section VI and Appendix C, that the averaging process produces a plausible, 
smoothly varying curve. 
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VI. ANALYSIS OF CURRENT DATA 



A. INTRODUCTION 

Various presentations were used to show different aspects of the 
current data, and to correlate this data with winds, tides, and ocean 
currents. A summary of results of the fourier analysis with mean current 
values and a table of contemporary current data are included in this 
chapter (Tables I and II) . 

The graphical presentations are contained in Section B of this 
chapter in the following order: 

(i) East-west current component versus wind and tide. 

(ii) North-south current component versus wind and tide. 

(iii) Total current magnitude versus wind and tide. 

(iv) Progressive vector diagram with the cumulative trans- 
port past the data point marked every five "fixes" 
(approximately 30 minutes) with a cross. 

Each set of graphs is prefaced with a short summary of sea and 
weather conditions for that particular period. 
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B. 



GRAPHICAL PRESENTATION OF DATA BY CRUISE 



1 . Cruise I 

Period: 10-11 May 1972 

Latitude: 36°46' N 
Longitude: 122°01'W 
Depth: 440 fathoms 
Ship's Speed: 5 knots 

Sea and Weather: Winds were from the WNW during this 
cruise, with a diurnal variation in wind speed caused partly 
by a land-sea breeze pattern. Velocities of up to 22 knots 
occurred in the late afternoon and virtual calm prevailed in 
the early morning hours. Seas varied from 0,5 to 2.0 feet 
with the wind velocity, while the swell was between 3 feet 
in the early morning and 8 feet in the early evening. Both 
were from the WNW. The sky was partly cloudy. Note 
that only the averaged current data is presented in this 
section. The unaveraged current data is presented in 
Appendix C. The start of the progressive vector diagram 
is labelled "S" and the end is labelled "E" . 
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Figure 5 a * Cruise I 
E-W Current Component, Tide 
Height, and hind, versus Time; 
Averaged Current Data 
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Figure 5c. Cruise I 
Total Current Magnitude, Tide 
Height, and Wind, versus Time; 
Averaged Current Data 
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Figure 5d. Cruise I 
Progressive Vector Diagram, 

N-S Component versus E-W 
Component; Averaged Current Data 
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2 . 



Cruise II 



Period: 30 - 31 May 1972 

Sea and Weather: Same pattern and magnitude as in Cruise I, 
but with fog and overcast skies in the morning, clearing by 
mid-afternoon . 

Note: All other information is the same as for Cruise l. 
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Figure 6d. Cruise II 
Progressive Vector Diagram, 

N-S Component versus E-W 
Component; Averaged Current Data 
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3. 



Cruise III 



Period: 19 - 20 June 1972 

Sea and Weather: Weak and variable winds with a maximum 
velocity of 8 knots were the rule, with a 1 foot sea and 3 
foot swell throughout the cruise. The sky was overcast and 
some drizzle fell in the morning. 
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Figure 7o. Cruise III 
N-S Current Component, Tide 
Height, and Wind, versus Time; 
AvoraRed Current Data 
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Figure 7d. Cruise III 
Progressive Vector Diagram, 

N-S Component versus E-W 
Component ; Averaged Current Data 
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4. 



Cruise IV 



Period: 22 - 23 June 1972 

Sea and Weather: Very similar to Cruise I with a maximum 
wind of 16 knots and overcast skies. Winds varied from 
west to southwest during this period. 
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5 . 



Cruise V 



Period: 10 - 11 July 1972 

Sea and Weather: Wind was fairly steady from the west to 
northwest at 5 to 18 knots with the same diurnal pattern as 
Cruise I. The skies were clear. Seas and swell were from 
the WNW with the same magnitude as Cruise I. 



52 




53 



Cruise 
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6 . 



Cruise VI 



Period: 20 - 21 July 1972 

Sea and Weather: Winds were predominantly southerly with 
velocities below 10 knots. There was a mixed swell, but 
it was mostly from the west and southwest with a southerly 
sea of 1 to 2 feet. It was overcast during much of this 
cruise . 
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Figure 10b. Cruiso VI 
N-S Current Component, Tido 
Height, and Wind, versus Time; 
Averaged Current Data. 
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Figuro 10c, Cruis® VI 
Total Current Magnitude, Tido 
Height, and Wind, vorsus Timoj 
Averaged Current Data 
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c. 



DISCUSSION 



1 . General 

This section deals with the mean currents and their relation 
to the wind and seasonal effects. The currents in Monterey Bay are 
weak and variable and have as yet defied accurate description. The 
current data gathered during this research and its relationship to oceanic 
currents, the winds, and seasonal effects may, together with other 
investigations, lead to a viable picture of the bay's circulation. 

The averages of the mean components for the six cruises 
are 11.6 cm/sec to the south and 3.6 cm/sec to the east. These yield 
a resultant of 12.1 cm/sec toward 163°T. No data are available which 
may be compared accurately to the present results. However, for purposes 
of rough comparison, Table I gives the mean for each cruise and currents 
derived from dynamic heights by Moomy (results yet to be published) . 

The points between which Moomy' s transports are computed are shown 
in Figure 4. The mean of Moomy's data is 15 cm/sec, and of the present 
data, 12 cm/sec. However, the standard deviations are relatively large 
and the close agreement probably is fortuitous. However, the generally 
southerly trend of all the mean directions of the time series data is in 
good agreement with Moomy's data. 

The "upwelling period," during which this current data was 
collected, has been described by Skogsberg (1936) as extending from 
February to September, and as having the same prevailing northwesterly 
winds as were experienced during this research. Lammers (1971) noted 
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TABLE 1 



COMPARISON OF RECENT OCEAN CURRENT MEASUREMENTS 



HOV/TON'S DATA 
MEAN SURFACE CURRENTS FROM 
24-HOUR TIME SERIES OF 
GEK MEASUREMENTS 
(AVERAGE OF 6 MEANS: 12 . 1/163°T) 


MOO MY’ S DATA 

TRANSPORT DERIVED FROM 
DYNAMIC HEIGHTS 






VEL/D1RC 






VEL/D1RC 


' DATE 


STATION 


CM/SEC/°T 


DATE 


STATION 


CM/SEC/ 0 








9 May 


3 & 5 


21 / SE 




Data 






3 & 4 


26 / SE 


10/11 May 


Point 


7.2 / 164 




4 & 5 


3 / ME 










2 & 4 


8/S 




(Near 






263 


27 / S 








16 May 


2 (Sc 3 


13 / N 


30/31 May 


Station IS) 17.1/175 


17 May 


2 & 3 


2 / S 


19/20 June 


Data 


4.1/198 


2 June 


2 & 3 


11 / S 


22/23 June 


Point 


11.3 / 13S 








10/11 July 




21.9 / 154 


6 July 


2 & 3 


22 / S 




Data 












Point 




19 July 


2 & 3 


7 /S 


20/21 July 




11.9 / 158 














25 July 


2 & 3 


14 / S 
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that the southern flow of the California Current exists at the seaward 
edge of Monterey Bay from mid-February to November, while the northerly 
Davidson Current persists during the balance of the year. Lammers also 
observed that during the upwelling period the mean velocity of the Cali- 
fornia Current in this region is from 1.0 to 6.7 nm/day with a mean of 
5 nm/day. This generally agrees in both direction and speed with the 
average of the mean currents from each cruise (5.9 nm/day). 

During this set of cruises the wind blew fairly steadily from 
the west northwest with the exception of Cruise VI, when the wind was 
from the south and southwest at ten knots or less. If the general 
southerly flow is due to the prevailing wind, this cruise would be ex- 
pected to show northerly flows. However, during Cruise VI the mean 
current maintained its mean southerly flow with only slightly more than 
the usual easterly component. In checking the progressive vector 
diagram for Cruise VI it was seen that almost due easterly current flow 
occurred during the daylight hours when the southerly wind was at its 
peak. A further study of weather charts for the period showed that the 
south wind was only a local phenomenon and that west northwest winds 
prevailed along the coast to the north. The local south wind did have 
some temporary effect on the surface currents, but the California 
Current maintained by the prevailing northwest winds along the coast 
persisted in its steady southerly flow across the mouth of the bay. 
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2 . Variability 

By inspection of the unaveraged Total Current Magnitude 
Graphs (Appendix C) it can be seen that the range of variability over 
24 hours is from 2 cm/sec to around 40 cm/sec. The Cruise V Total 
Current Magnitude graph (Figure 15c) is a good example of this variation 
and also shows two unusually high magnitudes of over 50 cm/sec. The 
spacings of the crosses on the Progressive Vector Diagram also give an 
indication of this current variability (each cross represents five current 
vectors, or approximately 30 minutes of time) as well as the changeability 
of the current direction. The Cruise I Progressive Vector Diagram 
(Figure 5d) shows the clockwise rotation (which may vary through more 
than 180 degrees) that is typical of all cruises except Cruise V. This 
directional oscillation is superimposed on a moderate mean southerly 
flow . 

In order to produce numerical estimates of the cyclic com- 
ponents of the currents and to determine their origins , the data was 
subjected to a fourier analysis. It was recognized that with so short 
a series (24.8 hours) the diurnal component would not be separated from 
the inertial component (20 hours at this latitude), and that if there were 
a substantial amount of noise it would have its usual effect in reducing 
the reliability of the coefficients. The N-S and E-W components 
separately were prepared by filling gaps, extrapolating to 24.8 hours 
total length, and by subtracting the mean. Fourier coefficients, mag- 
nitudes, and phase angles were then computed, using the Subroutine 
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FORIT from the IBM Scientific Subroutines Package (included in Appendix 
A). This program expresses the coefficients A and B as those in the 



equation 



m m 



f(t) = A + ^A cos m 60 t + S' B sin m 60 t (1) 

o m o m o 

where f(t) represents the time series being analyzed, m is the integer 

specifying the particular harmonic being calculated and XJ is the 

o 

fundamental angular frequency, which is related to the time duration of 
the series, T, by 



2 tt 



cj = _ 

o T 



The power present in the mth component is then 



( 2 ) 



P = 1/2 (A 2 + B 2 ) 
m m m 



and the magnitude of current speed in that component is 

, 2 2 1/2 
C = A + B ) 
m m m 



Equivalently, Equation (1) may be expressed as 



( 3 ) 



(4) 



f(t) = C cos(m CkJ t - (X ) 
m o m 



(5) 



where 



A =0 and OC 
o m 



sin Bm 
C 



(6) 



m 

Here OC is the phase angle, the number of degrees after the start of 
m 

each series at which component maximum occurs. 

The mean velocity and the phase differences between the diurnal 
component and the transit of the sun, and the semi-diurnal component 
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and the transit of the moon are contained in Table II, along with the 
results of the fourier analysis. The total power present in all the cal- 
culable frequencies is also tabulated in this table to facilitate a com- 
parison of this value with the amount of power present in the diurnal 
and semi-diurnal components. The conclusions to be drawn from this 
table are the following: 

(i) The amplitudes of the diurnal and semi-diurnal com- 
ponents are comparable with the mean velocity. 

(ii) The power in the diurnal plus semi-diurnal components 
is about three-fourths of the total variable power. 

(iii) Diurnal components are not consistently more important 
than the semi-diurnal and vice-versa. 

(iv) The variable effects are generally larger in the E-W 
components than in the N-S components, while for the 
mean velocities the reverse is true. 

(v) The diurnal and semi-diurnal terms in both the N-S 

and E-W portions of the current have a fairly consistent 
phase relationship with the meridional transit of the 
sun and moon respectively. 

The phase stability of the fourier components may be seen by 
comparing the phase angles relative to lunar and solar transit with the 
respective means. 

The variations from the mean of the semi-diurnal does not exceed 
+ 2.5 hours in either the E-W or the N-S component except for 
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,BLE II. FOURIER ANALYSIS RESULTS 
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Cruise III in the E-W and Cruises II and III in the N-S. Similarly, the 



variation from the mean of the diurnal phase does not exceed + 1.9 hours 
except on Cruises V and VI in the E-W and III and VI in the N-S. Only 
four out of the 24 variations from the mean exceed one-third of the asso- 
ciated period. Since we are dealing with a time series which is both 
non-stationary and short, such agreement may be regarded as good. 
However, the fact that there is not a constant relation between the 
amplitudes of these two components, and the fact that the phase angles 
occasionally vary widely from the mean leads to either of two inferences: 
(a) the currents are driven indirectly by the tides and are considerably 
modified by some unknown oceanic processes, or (b) inertial rotations 
are important. This makes one unwilling to say that the coupling of the 
tide and current is direct, but does indicate that a relationship exists. 
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VII. CONCLUSIONS 



Currents in Monterey Bay are quite variable in both velocity and 
direction. Therefore, one cannot describe the mean flow by means of a 
short-term measurement. 

All six cruises yielded mean southerly currents. The average of 
these means is a current directed 163°T at a speed of 12.1 cm/sec 
(5.9 nm/day) . This result supports Lammers' analysis that the California 
Current sweeps close to the mouth of Monterey Bay (and the data point) 
during the upwelling period. 

It has been established that there are substantial diurnal and 
semi-diurnal components which, with a few notable exceptions, have 
a consistent phase' relation with the sun and moon. There obviously 
must be at least an indirect coupling with the semi-diurnal tide. The 
diurnal component may be tidal or inertial or merely related to the passage 
of the sun . 
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APPENDIX A 



COMPUTER PROGRAMS 

1 . Reduction of GEK Data to Current Vectors 

2. Graphing of Ocean Current Data versus Winds and Tides 

3. Progressive Vector Diagram of Current Vectors 

4. Fourier Analysis of Ocean Current Components 
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APPENDIX B 



SAMPLE DATA FROM CRUISE I 
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CRUISE I 



10 - 11 MAY 1972 
WIND DATA 



TIME 
10 May 1972 



11 May 1972 



WIND SPEED WIND DIRECTION 



1230 


8 


NW 


1415 


11 


WNW 


1600 


12 


WNW 


1715 


16 


NW 


1800 


20 


WNW 


1955 


14 


NW 


2100 


10 


NW 


2300 


8 


NW 


2400 


8 


NW 


0200 


2 


VARIABLE 


0400 


4 


VARIABLE 


0800 


4 


VARIABLE 


0905 


6 


WNW 


1000 


7 


WNW 


1100 


8 


NW 


1200 


8 


WNW 


1300 


9 


WNW 
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CRUISE I 



10 - 11 MAY 1972 
TIDE DATA 

TIME TIDE I-IEIGIIT TIME TIDE HEIGHT 



10 May 1972 1200 


-0.90 


11 May 1972 2400 


LO 

CO 

O 

I 


1300 


-1.30 


0100 


-2.00 


1400 


-1.30 


0200 


-2 . 85 


1500 


-0.70 


0300 


-3.05 


1600 


+0.30 


0400 


-2.65 


1700 


+ 1 . 35 


0500 


-1.80 


1800 


+2.40 


0600 


-0.70 


1900 


+3.15 


0700 


+0.40 


2000 


+3.40 


0800 


+ 1 . 10 


2100 


+3.00 


0900 


+ 1.45 


2200 


+2.00 


1000 


+ 1.25 


2300 


+0.60 


1100 


+0.80 






1200 


+0.10 






1300 


-0.50 



/ 
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APPENDIX C 



UNAVERAGED CURRENT DATA PLOTTED VERSUS WINDS AND TIDES 
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Figure 12a. Cruise II 
E-W Current Component , Tido 
Height, and Wind, vs. Time ; 
Unavoraged Current Data 





97 



TIME (HRS) 




98 



Figure 13a. Cruiso III 
E~W.. Current Component, Tido 
Height, and Wind, versus Time 
Unavoragoct Current Data 
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Flgur© 13b. Cruise III 
N-S Current Component, Tide 
Height, and Wind, versus Tla© 
Unaveraged Currant Data 




100 



TIME (HRS) 




101 



TIME (HRS) 




102 



TIME (KRS) 




103 



TIME (HRS) 




104 




105 



V/NW NW WNW NW 




106 




107 



TIME (KSS) 




108 



TIME (HRS) 




109 



TIME (HRS) 




110 



BIBLIOGRAPHY 



1 . Barnes , C . A . , and Paquette , R . G . , Circulation near the Wash - 
ington Coast , Proceedings of the Eighth Pacific Science Congress, 
III, pp. 585-608, 1957. 

2 . Bowden , K. F . , Measurement of Wind Currents in the Sea by the 

Method of Towed Electrodes , Nature 171: 735, 1953. 

3 . Curtin , T . B . , Towed Electrodes in the Sea: Theory and Use , 

M. S. Thesis, Oregon State University, Corvallis, Oregon, 1970. 

4. Lammers , L. L. , A Study of Mean Monthly Thermal Conditions and 
Inferred Currents in Monterey Bay , M. S. Thesis, U. S. Naval 
Postgraduate School, Monterey, California, 1971. 

5. Longuet-Higgins , M. S., Stern, M. E., and Stommel, H., The 
Electridal Field Induced by Ocean Currents and Waves, with 
Applications to the Method of Towed Electrodes , M.I.T and Woods 
Hole Oceanographic Institution, Papers in Physical Oceanography 
and Meteorology, XIII (1), 1954. 

6 . McKay , D . A . , A Determination of Surface Currents in the Vicinity 

of the Monterey Submarine Canyon by the Electromagnetic Metho d , 

M. S. Thesis, U. S. Naval Postgraduate School, Monterey, 
California, 1970. 

7. Reid, J. L. Jr., A Comparison of Drogue and GEK Measurements 
in Deep Water , Limnology and Oceanography 3 , No. 2, 160-165, 
1958. 

8. Rodriguez, W. J., A Study of the Oceanographic Seasons at a 
Single Station in the Monterey Submarine Canyon 1951 - 1955, 

M. S. Thesis, U. S. Naval Postgraduate School, Monterey, 
California, 1970. 

9. Skogsberg, T., Hydrography of Monterey Bay, California, Thermal 
Conditions, 1929 - 1933 , Transactions of the American Philosophical 
Society, 29, 1936. 

10. Smith, T. D. , GEK Measurements of Surface Currents in Monterey 
Bay, 197.1 , M. S. Thesis, U. S. Naval Postgraduate School , 
Monterey, California, 1972. 



Ill 



11 . 



Stoddard , H . S . , Feasibility Study on the Utilization of Parachute 
Drogues and Shore-Based Radar to Investigate Surface Circulation 
in Monterey Bay , M. S. Thesis, U. S. Naval Postgraduate School, 
Monterey, California, 19 71. 

12. Sysoyev, N. N., and Volkev, V. S., A Guide to an Electromagnetic 
Method of Measuring the Velocity of Sea Currents from a Ship Under 
Way, Trudy Institituta Okeanologlii, Vol. XXIV: 173-199, 1957. 

13. U. S. Department of Commerce Coast and Geodetic Survey, Tide 
Tables - West Coast North and South America , U . S . Government 
Printing Office, Washington, D. C., 1971. 

14. U. S. Department of Commerce Coast and Geodetic Survey, Tidal 
Current Tables - Pacific Coast of North America and Asia , U. S. 
Government Printing Office, Washington, D. C., 1971. 

15. U. S. Naval Oceanographic Office, Instruction Manual for Obtaining 
Oceanographic Data , Pub. No. 607 , U. S. Government Printing 
Office, Washington, D. C., 1970. 

1 6 . von Arx , W . S . , An Electromagnetic Method for Measuring the 
Velocities of Ocean Currents from a Ship Underway , M.I.T. and 
Woods Hole Oceanographic Institution, Papers in Physical Ocean- 
ography and Meteorology, XI (3), 1950. 

17. Wyllie, J. G., CALCOFI ATLAS No. 4, Geostrophic Flow of the 
California Current at the Surface and at 200 Meters , State of 
California Marine Research Committee, 1966. 



112 



INITIAL DISTRIBUTION LIST 

No. Copies 



1. Defense Documentation Center 2 

Cameron Station 

Alexandria, Virginia 22314 

2. Library, Code 0212 2 

Naval Postgraduate School 

Monterey, California 93940 

3 . Department of Oceanography 3 

Naval Postgraduate School 

Monterey, California 93 940 

4. Dr. R. G. Paquette 5 



Department of Oceanography 
Naval Postgraduate School 
Monterey, California 93940 

5. Dr. J. J. von Schwind 1 

Department of Oceanography 

Naval Postgraduate School 
Monterey, California 93940 

6. Dr. W. C. Thompson 1 

Department of Oceanography 

Naval Postgraduate School 
Monterey, California 93940 

7 . Dr . D . F . Leipper 1 

Department of Oceanography 

Naval Postgraduate School 
Monterey, California 93 940 

8. Dr. E. B. Thornton ' 1 

Department of Oceanography 

Naval Postgraduate School 
Monterey, California 93940 

9. Dr. N. E. J. Boston 1 

Department of Oceanography 

Naval Postgraduate School 
Monterey, California 93940 



113 



No. Copies 



10. Dr. W. W. Denner 1 

Department of Oceanography 

Naval Postgraduate School 
Monterey, California 93940 

11. Chairman, Association of the Monterey Bay 2 

Area Governments (AMBAG) 

Monterey, California 93940 

12. Oceanographer of the Navy 1 

The Madison Building 

732 N. Washington Street 
Alexandria, Virginia 22314 

13. Lieutenant H. M. Howton, USN 3 

149 Ferry Road, N.E. 

Fort Walton Beach, Florida 32548 



114 



ral Postgraduate School 
nterey, California 93940 



essra»9*3r*nQca 



Securit y Cla ssification 

""document control data - r & d 

(Sfcurj ty classification of title, body of abstract end indeylng annotation must be entered when the overall report is classified) 
cTnaTING activity (Corporate author) 2*. REPORT SECURITY CLASSIFICATION 



Unclassified 



26. GROUP 



ORT TITLE 

tudy of Time Variability of Surface Currents at a Point in Monterey Bay 



SCRIPTI VE NOTES (Type of report nnd 9 Indus ive dates) 

ster's Thesis; December 1972 



TKORISI (First nsmo, middle Initial, last nemo) 



rry M. Howton 



ORT DATE 

;ember 1972 



7*. TOTAL NO. OF PAGES 



116 



76. NO. OF REFS 

17 



)N T R A C T OR GRANT NO. 



*OJEC T NO. 



ORIGINATOR’S REPORT NUWBER(S) 



Pb. OTHER REPORT NO<S> (A ny other numborn that may be err/^ned 
thle report) 



STRIBUTION STATEMENT 



moved for public release; distribution unlimited 



IMPLEMENT ARY NOTES 



12. SPONSORING MILITARY ACTIVITY 

Naval Postgraduate School 
Monterey, California 93940 



3STRACT 



The geomagnetic electrokinetograph (GEK) was used to measure surface currents 
ar the center of Monterey Bay during six separate 24 -hour periods from May through 
y, 1972. An average of 244 current vectors were derived for each cruise. The mean 
rents from these cruises are all southerly and ranged from 4.1 cm/sec to 2 0.4 cm/se 
3 average of these mean currents is 12.1 cm/sec toward 163°T. These values were 
npared with individual currents derived from dynamic topographies from the some 
•iod. Diurnal and semi-diurnal variations of the current were studied after subjecting 
data to a fourier analysis. It was concluded that there must be at least an indirect 
jpling of the ocean currents with the semi-diurnal tide at the data point. The diurnal 
mponent also is important; it may be tidal or inertial, or merely related to the passage 
the sun . 



) Ft 
9 1 N 



FORM 



ov # J473 

0101 -807-C8I 1 



(PAGE I) 



1 3 5 



Security Climii flea lion 



A- 3 1 408 



Security Classification 




o 1 o 1 • e o 7 - 6 e 2 1 Security CU»*ificetion a-3M09 



The. 

H82> 

c.l 



Thes i s 

H8265 

c.l 



■t r 

- v 



soi: 



Howton 

A study of time vpr- 
iability of surface 
currents at a jn 



OCTMcrtiterey Bay. 



ch ‘ j A I: '(it) 
6 JAS 67 



2 27 5 n 

2 5 7 0 0 

S I 0 3 3 4 



Thes i s 

H826t 

c.l 



r ■* 

•< 



301 



Howton 

A study of time var- 
iability of surface 
currents at a point in 
Monterey Bay. 






I -1 % . <»ff • 



II Mill I Nil 

3 2/G8 00? 00/02 1 
chj> 1 1 >H0‘ uonAn» 




